home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zhpgst.z / zhpgst
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZHHHHPPPPGGGGSSSSTTTT((((3333FFFF))))                                                          ZZZZHHHHPPPPGGGGSSSSTTTT((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZHPGST - reduce a complex Hermitian-definite generalized eigenproblem to
  10.      standard form, using packed storage
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE ZHPGST( ITYPE, UPLO, N, AP, BP, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, ITYPE, N
  18.  
  19.          COMPLEX*16     AP( * ), BP( * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      ZHPGST reduces a complex Hermitian-definite generalized eigenproblem to
  23.      standard form, using packed storage.
  24.  
  25.      If ITYPE = 1, the problem is A*x = lambda*B*x,
  26.      and A is overwritten by inv(U**H)*A*inv(U) or inv(L)*A*inv(L**H)
  27.  
  28.      If ITYPE = 2 or 3, the problem is A*B*x = lambda*x or
  29.      B*A*x = lambda*x, and A is overwritten by U*A*U**H or L**H*A*L.
  30.  
  31.      B must have been previously factorized as U**H*U or L*L**H by ZPPTRF.
  32.  
  33.  
  34. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  35.      ITYPE   (input) INTEGER
  36.              = 1: compute inv(U**H)*A*inv(U) or inv(L)*A*inv(L**H);
  37.              = 2 or 3: compute U*A*U**H or L**H*A*L.
  38.  
  39.      UPLO    (input) CHARACTER
  40.              = 'U':  Upper triangle of A is stored and B is factored as
  41.              U**H*U; = 'L':  Lower triangle of A is stored and B is factored
  42.              as L*L**H.
  43.  
  44.      N       (input) INTEGER
  45.              The order of the matrices A and B.  N >= 0.
  46.  
  47.      AP      (input/output) COMPLEX*16 array, dimension (N*(N+1)/2)
  48.              On entry, the upper or lower triangle of the Hermitian matrix A,
  49.              packed columnwise in a linear array.  The j-th column of A is
  50.              stored in the array AP as follows:  if UPLO = 'U', AP(i + (j-
  51.              1)*j/2) = A(i,j) for 1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-
  52.              j)/2) = A(i,j) for j<=i<=n.
  53.  
  54.              On exit, if INFO = 0, the transformed matrix, stored in the same
  55.              format as A.
  56.  
  57.      BP      (input) COMPLEX*16 array, dimension (N*(N+1)/2)
  58.              The triangular factor from the Cholesky factorization of B,
  59.              stored in the same format as A, as returned by ZPPTRF.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZHHHHPPPPGGGGSSSSTTTT((((3333FFFF))))                                                          ZZZZHHHHPPPPGGGGSSSSTTTT((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      INFO    (output) INTEGER
  75.              = 0:  successful exit
  76.              < 0:  if INFO = -i, the i-th argument had an illegal value
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.